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ABSTRACT 



False- positive and false negative decisions are the fundamental 
errors comnitted with a mastery test; yet the estirnation of the - 
likelihood of coffsiiltting these errors has not been investigated. 
Accordingly, two methods of estimating the likelihood of committing 
these errors are described and then investigated using Monte Carlo 
techniques. Conditions for obtaining accurate estimates are noted. 
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ESTIMATING THE LIKELIHOOD 0^ FALSE-POSITIVE AND 
■ o ■. FALSE-NEGATIVE DECISIONS IN MASTERY TESTING: 

AN EMPIRICAL BAYZS APPROACH ' - 

1. ■ ' Introduction - 

Typically a mastery test is designed to sort k examinees into one of two 
mutually exclusive groups. For example, in a program of Individually Rres-. 
cribed Instruction a student's progress through^each level of a program of 
study is governed by his performance on a test dealing with individual be- 

o , 

, havioral objectives. The. purpose of a test in such situations is to make 
a mastery>(nonmastery decision for each of k examinees* If a mastery deci- 
sion is made for a particular examinee. then he is advanced to the next level 
of instruction.. If, however » a nonmastery decision is made he w,ill be given 
remedial work. 

A model of mastery testing which is frequently adopted may be described 
as follows: A pool or domain of dichotomously scored test items, having 
mixed item difficulty* is constructed in relation to a particular course of. 
instruction. The item pool may exist de;: facto or it may be a convenient 
conceptual i^zationi The item form no^" jvely and others (1973) repre- 

sents such a conceptualization. Let ue the percent correct domain 

th ^ 

or "true" score of ^tlje i examinee; A^. represents the percent of .Items that 
the ith examinee would answer correctly if he were to respond to every item 
in the item pool at a given occasion in time. With respect to the domain of 
items an examinee is said to have attained mastery if x. >^ and nonmastery 
if A. < A where' A is a known constant with a value between zero and one. 
The problem is to make a mastery/nonmastery decision for a given examifiee' based 

/ ^ 

o ■ . . . . 
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on 'his responses to n items nandomly selected from the item domain. A mastery 
decision is made if an examinee answers n^ or more items correctly where ^ 
0 1 ^^o 5 Note that corresponds to a'concept of mastery and n^ to a 
mastery score or index (Harris, 1974). We further observe that the model of 
mastery testing just described Is equivalent to the ranking and selection pro^ 
blem of partitioning k populations (examinees) with respect to a standard. 

This model provides^ a reasonable description of mastery testing and is 
consistent with definitions of mastery or criterion-referenced tests (Glaser 
and Nitko, 1971; Harris, 1974). (See also Hambleton and Novick, 1973; Fhah'er, 
1974; Novick and Lewis, 1974; Huynh, 1976; Wilcox, 1976). 

A false-j^oslTive error occurs when the exajniner estimates ah examinee's/ 
true score to be above the criterion level Xq when in fact it is not, A( 
false-negritive error occurs when x. is estimated to be below Xq when the reverse 
is true. False-positive jand false-negative errors are the two errors that can 
be made in a two-valued classification, yet the estimation of the probability 
of committing these types of errors in connection with mastery testing has 
been virtually ignored. Instead attention has ^een given to measures of stabi- 
lity such as the proportion of agreement (Hambleton and Novick, 1973) which 
estimates the probability of randomly selecting an examinee and classifying 
him the same way based on two administrations of the same test. Certainly it 
is desirable to have a test with a high degree of stability. However, it/may 
bff that such attest, is consistently inaccurate. 

Let a and e equal the probability of committing a false-positive and false- 
negative decision, respectively, for an examinee chosen at random from some 
population of potential examinees. -Observe that the values of both a and e are 



a function of the number of items^on the test as well as the instructional . 
history of the examinees. Consequently, knowing" a and g provides a meaningful 
characterization of the entire teaching-testin^complex. The purpose of this - 
paper is to examine the problem of estimating a and 6 based on student response 
data. 

The binomial error model gives a reasonable approximation to the observed 
score distributions- on tests (Lord, 1965, p. 253), but the compound binomial 
may be more realistic (Lord, 1965, Section 6; Lord- and Novick, 1968, Chapter 
23) and hence may gTve more accurate results. Accordingly, two methods of 
estimating a end 8 are described and the accuracy of these statistics are 
examined under both the binomial and compound binomial error modf- . Sections 
3; and 4 derive estimates of a and e assumlntj that. the distribut/ion of true 
scores belongs 'to a particular parametric family. - Section 5 examines the 
accuracy of these estimation procedures using Monte Carlo techniques. 
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'Z: Mathematical Statement of the Problem - 

As indicated earlier, we let \^ denote the proportion correct "true", score 
of an examinee taking an n item test. We regard x. (i=l,...,k} as a sample 
from a prior distribution, say g(x), where'O < x <. 1. Let h(x | x.) be the 

distribution of observed scores for a gjven true score x.. Assuming that g(A) 

. . ■ ■ t 

Is an integnable function and since h(x | x.) is discrete we havej 

V (2.1) ,a == r /o° h(x Hx) g(x) dx. , ■ . j^--^^^^^^ ' 




; and - 

no-T , - ' 

(2.2) 6 = S -/J h(x I X-) g(x) dx 

x=0 



If we knew g(x) and if we assume h(x^ I x) is binorularwe would also know 
a and e. However* g(?^) is usually unknown. The approach taken here is to use 
empirical Bayes procedures to estimate g{X). When 4i(x | x) is assumed to be 
compound binomial, it too must be estimated. We are particularly interested 
in the accuracy of point estimates of a and e for relatively small values of 
k and n. Me emphasize that the results given below do not reflect directly 
the accuracy of our estimate of g(x). In fact we are 'only concerned with an 
accurate estimate of g(x) in so far as It improve^ bur estimate of a and 6. 
It may be» for example, that a relatively ppor estimate of g(x) will yield 
a reasonably accurate estimate of the- frequency of occurrence of both false- 
positive and false-negative decisions. ^ 

As a measure of the accuracy of any statistic « wh^ is used to ^esti- 
mate a we use the expected value of the square of the difference of a and 5 
over the joint distribution of x and x\ That is, we use 

(2.3) w^ ^ S fl (a-S)2 h(x I X) g(x) dx. 

which corresponds to the average risk used "^n empirical Bayes methods (see, 
e.g., Maritz, 1970, p. 3). When estimating i.. with 6 we use w^ which is 
defined by replacing (a-a) with (3-0) in (2.3). 

3. Estimation of a and BcassCiming a beta prior 

The es'timation of the prior distribution g(A) is a most difficult problem 
for which no general solution exists. It behooves us, therefore, to consider 
the estimation of g(x) under a variety of conditions. We begin with perhaps 
the simplest, but also the most severe restriction on the family of prior 



olstn'butions, namely, that g(x) is an incomplete beta distribution with 
t)aranieters r and s. That is, /x 

/.X r (r + s) .r-1 /, ,*s-l 

where r is the usual gairnia function. This assumption is restrictive because 
there is little if any doubt that g(A) does not belong to the family of beta 
priors; Yet th^re are several reasons for considering this case. First, the 
beta density is the natural conjugate prior of the binom,ial kernal 'h(x | A) 
(Raiffa and Schl.aifer, 1961). Second, we have indentifi ability, i.e., there 
exists a unique g such that 

\ (3.1) f(x) = /J h(x U) g(^) ' \ 

where f(x) is the marginal distribution of observed scores (Mciritz, 1970. 
chapter 2). Third, there is evidence that a reasonable though not entirely 
satisfactory approximation of the true score distribution canjbe obtained with 
a two parameter beta prior (Keats and Lord, 1962). Finally, and .perhaps most 
importantly of all, the results of estimation procedures assuming a particular 
parametric form for the prior may be used as a bench mark for j;dging alter- 
nate estimation techniques. 

Let X.. (i=l,...,k; j^l,...,n) denote the j*^ observation on the i*^ , 

. - \ ' ' ''' ' ' 

examinee". By estimating the parameters r am s of the beta prior based on the 

sample x.., we obtain an estimate ^of g whi|:hUn turn yields an estimate of 

both a and b- We begin by describing a method of estimating r and s which 

assumes that the binomial error model (Lord and Novick, 1968, ^chapter 23) holds, 

h(x I X) =(j;)AMl-x)""' ' 



X 



Let Mj-^-j denote that t factorial moment of the marginal distribution of 
observed scores, i*e. , 

Let y. represent t]iG t moment of the true score distribution g(x). Then 



(3.2) = 



(n - t) !M 



t nV 

(Lord and Novick, 1968» expression 23.8.4). We ob.tain unbiased estimates of 
Mp-j and Mj-g-j with 

(3.3a) Mr.T = l/k T.^ n. ' 

(3.3b)" Mj-23 = l/k xlj,^ (x? - x.) 

where x^ = s"^^ x^j. From (3.2) we have that 

(3.4a) ■ y = Mj-'^-j /n . " ' _ 

> 

(3.4b) > {ig ^ '^[2] ^ - I)) • ' 

\ . ■ / ^ 

are unbiasedl estimates of and Thus^, the mean, say u and the variance, 

PI 

say a , of the prior distribution may be estimated as 



(3«5a) i y * y 



1 



(3.5b) = {Ig-^^ " 

Note that it is possible to have Mj-gj ='0 and M^j > 0 resulting in a negative 
estimate" of a^. When this. occurs we estimate each a. as , 

• k n 



e'\ x../(kn) 
i=l j=l 



A . \ ^ 

For this special condition, if X >.„Xo we estimate a to be zero and B to be 

^,,0 ^x^^"^^ (14)""'^. Correspondingly, 1f X < Xq, wg estimate 6 to be zero 
and o to be eJ^^^ (Hj {if (1-a)"-,\ 

In terms.of r and s . ' * 

(3.6a) . - r 



r + s 



(3,6b|) -2 _ , 



(r + s)2 (r + s + 1) 



(3.7a) r = 0 ~ u) _ ^ 



G 



(Johnson and Kotz, 1970). Solving (3.6) fdr r and s gives ' 

. (3.7b) s = ii-X^-Jil + p -1 

Substituting In (3.7) u and a for u and a , respectively, yields an estimate \. 

of r and s, say f and s. The estimates of r and s may then be used to 

estimate a and g as , ' ^ 

. ' (3.8a) a. = ? f,lo (n) x>< (3.x)n-x ilL±Jl ,r-l(^^,jS-l 

x=n^ r(r) r(s) 

(3.8b) = ^'^ (5) (1-x)"-'^ A^-''(l-x)^-^ dx 

' x=0 ^° ^ r(.r) r(s) 

As indicated earlier, ttie binomial error model may not be completely satis- 
factory in an item sampling model • As suggested by Lord anid Noyick (1968, 
chapter 23) -^e use a two-term approximation to the compound binomial , viz.. , 



(3.9) h(x I X)' = (x) + dx(l-x) C(x) 



where 



Lord (1965) notes that (3. &) is a close approximation to a f*-e^^uency distribu- 
tion for most cases of interest. Difficulties could arise if d were too large; 
we avoid these difficulties by assuming 0 <.d < 4.0. For all 16^distributions 
reported by Lord, the values of d were in this„range. , (See Lord, 1965^ p. 264). 

Under the more general compound-binomial we are still able to estimate a , 
and e.-- As shown by Lord (1965, p. 265) the mean and variance of the distribu- 
^tion of true scores for the two-term approximation to the compound binomial 
error model are given Iby: 
\ (3.10) Pi = Mp-j/n 

\ (3.11) , 0^ = al - (n-2d) p q 

■ \ ?. i ■ ■ ■ - 

where c is the varic^nce of the distribution of observed scores, p = Mp-j/n and 



X 



q = 1-p. The parameter S i^ given by 

V n^ (n - i'^- a^- ^ . 

(\3.12) d = ^ ^ 

\ 2[mj^ (n - M^) - - n a J ^ 

where is the variance of the item difficulties, d may be 'estimated u^sing 
standard\i tern analysis techniques'. Hence, i.he paramstsrs r and s of the ^eta 



Substituting (3.9) into {3.8a) s the estirnate, of a under the" compound ' 
binomial error model is 



, A ^ , 



., ^3.13a) ./^i{x V-^)^^fc^ x""^(l-^)^'' c!^ 

x=no r(r) r{s) 

. • • , ■■ ■ ■ ■ V •/ 

Correspondingly* we estimate e to be / \ | 

(3.13b) S,,„ /J h (X I .) liLtii-' dx 

KO} Ao . rtr.).r(s) ; 



4. Estimation of a and B using an invers^ sine transformation. 

In the V^vious section a procedure for estimating ct and 8 was deiscribed 
which is. contingent upon estimating the parameters, r and s of an assumed beta 
prior. One difficulty with this estimation procedure is that the sta^tistics r 
and s no doubt- lack the desirable properties of^unbf^^edness, maximum likeli- 
hood, arid efficiency. Consequently, one might expect estimates of rl arid s to 
be poor i\or relatively smaTl samples'. Since one^ would, hope that accurate esti- • 
mates of r\and s would yield accurate estimates of/t and' e it may 'fee helpful to 
search for more accura^e^estimates of r and s even though improveipent in our 
estimates of i^and s promises to be a most diffiTGUl^t task. For, ex^imple, even 
in the simpler Lore conventional case in which the , sampled' values are known, 
maximum likelihood estimates ot r and s are obtained iteratively.r We propose, 
therefofe.lto investigate the u«e of an inverse sine transform wnicli converts 
a binomial random variable into an approximately normally distributed random • 
variable with known variance, the variance being independent of the value of 

— — ' • ■ N * ' 

this is often called a variance „sVabi 11 zing transformation. The advantagKof this 



approach is that the est1 characterizing th- r 

distribution can be expected to be more accurate rsl a l^ive to"tli^eta-L nomial ■ 
model described ab'ove if the transformation used does indeed yiel^a normally 
distributed random variable. The disadvantage of this approach Isvthat the 
transformed random variable is asymptotically normal and thus any estimation , 
/ ./procedure usihg small samplS may be poor. In additionv the rate of convergence 
' to normality is a function of the unknown parameter x^. The crucial question is 
of course, whether this approach reduces the values of and as defined l?y 
(2.3) above. ^ 

(4.1) y^ =■ (4n + 2)^ (sin"^ (/_^) sin"'' (/""^Trr)) 

where, as before, » s x.. is'the observed score of the i examioee. 

d=l . c 

, The transformation (4-1 ) is Suggested by Freeman and Tukey (1950) where y^ 

' ' ' ■ ' \ ■ ' * ■ • 

givert >j is approximately normal ly distributed with mean 

• ■ / >■ ■ . \..-^ ■ . ■ . 

. t ;= (4n + Zr sin"' (/x^) i 

and variance one. As i ft the previous section we let )) and represent the 
mean and variance of the prior distribution. Herer^owever, the natural con- 

'. jugate prior has a normal distribution (Raiffa and Schlaifer, 1961). Moreover, 

■ • , V ''■ ■ ■- ' ■ ., • 

i ■ the marginal distribution of observed score is also normally distributed with 

. ■ - --"^^ ■ 2 ■ ' 

meari^K and* variance a + 1. It follows that 

: ^E^ E{y | r) = E(t) . ■ ^ . ■ 

v.:vf ; ■ " i E^ E(y^ - 1 I t) = ' .' 

K'r we may estimate y = E(t) and a-.? E(t ) - E (t) with 
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■ 1 k \ 

(4.2a) G = 4- E • . 

i=l ' 



It is. known (see^ for example, Hogg. and Craig, ,1970^ 210) that the joint 
probability density function of y and x is biyariate' no rmal ^w ith common mean 
• respective variances 1 + and o^, and correlation a /.A + o^. Applying the/ 
i;^thod,bf moments,. as Was done in the previdus section, the estimates of a - 
and e are ^ ' 

: . (4.3a) -z'' fy f^^ exp i^^^ 

•2'a ■ J- ' 1-' * 



where 



and 



A' n (4n+bs1n",\(/x") . 

M ■ ' 

Again we have the difficulty that may be negative. In this case we set 

/ * n 1 ' * . ' ' * 

Op = 0/ B« = (5) x'^(l-x)""^ when X >. XqI and when X < x^ we set e^'O, 

I ,x=o- ;-/ 

ctp = 2 (5) X (1-a)' as was done in the previous section, 

. * It may be helpful, to indicate how (4.3a) and (4.3b) can be evaluated with 

I.. . ' " 

^ existing computer subroutines. For convenience we write the bivariate 

■ •' ' . ,'11 , 
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distribution of y and A as h(x | X) g (x). Observe that 



(4.4) 1 X) g(A) d\ dy+/l^ h(x | xj g (x) dx dy 



where f(y) is the normally distributed, marginal distribution of the observed 

scores y with mean p and variance 1 + o^. The first integral on. the left hand 

■ ^ • • • ■ ' . ' ^ . ... ' _ „ • ■•• 

side of (4.4) can be evaluated with the IMSL subgouj5LneJli^^ 



the random variables y and-X are trans;fonned so as to have common mj?an-^^r^^^ 

and variance one. The right hand side of (4.4) canjfee^evaltiate^^ the FORJRAN 

subroutine ERFC J thus, we have the value-g^^T^e^ statistic may be evaluated 

■ \ "^"'^"^ ■ ,/ ^" "' . ■ . ^ ^ \ 

in a similar manner. ' ^ , 

The transformation (4.1) is claimed by Mpsteller and Youtz (1961) as well 
as ^josteller and Tukey (1968) to be the best existing angul air transformation 
for the'' binomial distribution. As indicated above, however, the compou^id 
binomial may be a more appropriate probability distribution for descrtbingVthe 
observed frequency of test scores. To be conservative we introduce an inverse 
sfne transformation for the compound binomial. The desirability of using this 
transformation will be discussed in section 5 below. ^ 

As^ shown by Lord (19.65, p. 266) the mean and variance of the ttyo term 
^approximation to the compound binomial distribution given by (3.12) may be 
written as nx and (n-2d) x (1-x). It follows that 

V 7 -.-^ x. - 



\ 



(4.5) , >-2n_ sin""" (vOT) . 

I- ■ .. . . ^ . : 

is asymptotically normal with mean ^" _ and variance one (Rao, 1973, 
^ ^ n-2d 

, ^ . ':. ■ 15 ■ ■ ' 



Section 6g). After estimating d via (3.11), one may use transformation C4.5) 



in place of (4.1), then estimate v and with (4.2a) and (4.2b); respectively, 
anr' f> a and e with (4.3a) . > i (4.3b). When a and 6 are esti- 

mated using (4.5) for a given value of d, we denote the estimates by Cg^^j 

and'62((ij» respectively. . ' 



5. Method and Results of Monte Carlo Experiments 

The true i'tore for each of the k examinees (the value of x^, i|=l,. ..,k) 

was generated according to a beta distribute u with parameters r and s. The 

priors used thcluded i-shaoed, U-shaped, sy, letric and skewed disjributions.. 

Some these distrtbut is (e.g. U-shapfed) • ?robably unrealistic in 3nns 

of rrsc- ary testing. The ; were included, howr so as to obtain; more genera! 

resi 9. Onc(? x. was determined the observed jre was generated accord- ng 

to the two-term approximation to the compound binomial given by (3.9) for 

d=0.0, 2.5, 4.0. Expression (3»9) was evaluated by using the relationship • 

I (a, n-a+1) = 2 \^ (l-x)""^ in cotijunction with IBM's SSpfll^) subroutine 
^ x=a ■ ■ ■ ; r~ . V <• , ' [ 

BDTR whare 1^^ denotes the incomplete/bdta function ratio (see dphnson af\d Kotz, 

19'0, chapter 24). Over 200'Monte Carlo studies were made for/each of the 

/ ' - / • ■ • 

es'^fp' ors oil, 3i, 02 antt 82- . , ,• I 

Initially we set Xg=0.7 andji^=x^n. For eacK^prior disti|ibution used, 

w ur.: w. were estimated jy first/using the exact value of d ^ind then by 

a b ' . / • ■ \ - 

setting d arbitrarily equal to z^ro. The values of k and n Were (k, n) = 
(10, 10)," (10, 20), (10, 30), (20. ,10), (20, 20), "(30, 10). jAll estimates • 

of w and Wu were based on 500 iterations. For simplicity v>/e discuss the 

, a D , ■ , . / ■ . • i 

results in terms of w^. No additioral insights were found when examining Wj^. 



" Regardless of the true score distribution ustd, the value of d had 
negligible effect on the value of when (3.8a) was used to estimate a. 
This result is illustr-'ad in Table I for the special case (r, s) ' (9, 2) 
and (r, s) = (3, o). = .7. Moreover, using the exact value of d in 
(3.13a) generally had little effect on lowering, as^ demonstrated in Table 
II. One exception to this finding occured for k=n=10 arid r=s=3. For d=0. ~ 
w was est^mated%o be .093 using {3.8a). For d=4.0, w dropped to .079. 
In gener , however, varytig the value-of d affected w^'only at the third 
decimal This result also held when the more general (3.13a) was used 

which incit-ucrats" the, two- term approximation to the compound ainomial. This 
finding was net surprising since there- appe?rs to be negligible change in the 
observed Sxcore ^'.c true score distributions when' the value of d is altered 

(Lor-. 1-^% . " . 

As a,, again it was found that altering d had little effect the 
value of ■ However, setting d arbitratiiy equal to zero and using tra.r.s- 
.formatior . : tended tb give better res.Tts (lower values of Wg) as opposed 
to usinc tne e^ ct value of d and transfot-ation (4.5). The apparent reason 
for this c i is that (4.1) converges more rapidly to normality than does 

(4.5). ^ 

Table I a -sil as Table II suggests that reasonably accurate estimates 
of a can bt m sained particularly if n is greater than or equalto 30. To 
better assess tne accuracy of ai and' aj we^resent Table III which gives the 
values of or both aj and where r=9, s=2, x^=.,5(.l).8. From T^ble III 
we see that c and ^2 a^e very accurate for A^=.5 but that this accuracKdi^mi 
nr-ihes c«- ^fe^ably as approaches .8. The difficulty is that both «i an 



tend to underestimate o. The results Indicate that the amount by which a 
is underestimated increases as a gets large. In Table III, for example, the 
actual value of a is 13.005 when k=n=10 and x^=.5. For x^=.8, a=0.157. 
V Note that for x =.8 the most effective method of lowering is to 

increase n (the, number of items) as opposed to increasing k (the number of 

■ ■ ■ ■ . ■ - ..' ■« ■. • ■ - ' " ■ ' " . ■■■ 

individuals). In general , but not always, increasing n will decrease the value 

of o. Consequently, /we lower the value of w, by_ Increasing n primari ly because 
We obtain more Accurate estimates Of a when a is small. Increasing k with n 
fixed also lowered w^ but at a much slower rate. 

We also observed that heiti^er of the statistics Si or Sa dominated the 
Other, i.e., had consistently lower values for w^. Consequently, based purely 
on statistical considerations, 1t*is impossible to recommend one method of 
estimation rather than the other. However, we see that ai dominated 02 for 
.7^ ..8 particularly for = .8. The reason „is that S2 = 0 occured 
more frequently due to oegaiive estimates of the variahc^ of the prior. Since 

the values of w were particularly large for x^ = .8, it would' seem best to 

■ .a - . ■•■ - , ° •■ 

use ai-. On the otherhand, to ensure, accurate estimates of a, it would 'seem . 

prudent to have n equal to at least 30 and preferably largek In this 

case, evaluating Si might be more dif^ficult cpmputationally. By having n large 
however, accurate es.tim|ites may still be possible with Sa* 



\ 



TABLE I 



Values of using (3".8a), Xq-.? 





10, 10 


10, 20 


MO. 3| 


20, 10 


20, 20 


30, 11 


- 






> 

s=2 








0.0 


.043 


.035 


. ;;p31 N-^ 


.036 ^ 


.028 


.033 


2.5 


,. .041' 


.035 


.p31 


.036__/. 


.029 


.034 


4.0 


.041 


.035 


.031 


.036 


.030 


.034 








r=3, s=3 








0.0 


.093 


.059 


.045 




.052 


.087 


2.5 


.083 


.057 


.044 


.oar- 


.051 


^081 


4.0 


.079 


.056 


.043 


.077 


.051 


, .077 
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. TABLE II ' " 

Values of using the exact value \. 

' ' of d in (3J3a), \o^J 

10.20 ^ 10 . 30 20, 10 20 . 20 30, 10 
r=9, s=2 

. * ■ ' ' \ 

.034 . 029 .037 \029 ° .035 

.034 .029 .036 .030 .035 

.052 .038 .082* .051 .082 

.051 .03g .080 .060 y .077 



( 
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, TABLE - il 

Values of w,, r=^, s=2, cl=0 
a 



10» 10 10, 20 10, 30 20. 10 . 20. 20 
Results using aj 



30, 10 



.017 


^XOII 


.008 .016 


.020 


J 

.014? 


.012 .018 




.035 


.031^ .036 


.106 


.082 


.070 .094 






Results using 02 


.006 


.004 


.004 • .004 


.019 


.016 


.014 .020 


/056 


. .035 


.041 ..058 


. 1 27 


.100 


.085 .127 

■ ■ , ■ a 
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